Statistical mechanics of 2D turbulence with a prior vorticity distribution 
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We adapt the formalism of the statistical theory of 2D turbulence in the case where the Casimir 
constraints are replaced by the specification of a prior vorticity distribution. A phenomenologi- 
cal relaxation equation is obtained for the evolution of the coarse-grained vorticity. This equation 
monotonically increases a generalized entropic functional (determined by the prior) while conserving 
circulation and energy. It can be used as a thermodynamical parametrization of forced 2D turbu- 
lence, or as a numerical algorithm to construct (i) arbitrary statistical equilibrium states in the sense 
of Ellis-Haven- Turkington (ii) particular statistical equilibrium states in the sense of Miller-Robert- 
Sommeria (iii) arbitrary stationary solutions of the 2D Euler equation that are formally nonlinearly 
dynamically stable according to the Ellis-Haven- Turkington stability criterion refining the Arnold 
theorems. 

PACS numbers: 47.10.A-,47.15.ki 



I. INTRODUCTION 

Two-dimensional incompressible and inviscid flows are 
described by the 2D Euler equations 



'dt 



u- = 0, 



-zxVV', (1) 



where u is the vorticity and ijj the streamfunction. The 
2D Euler equations are known to develop a complicated 
mixing process which ultimately leads to the emergence 
of large-scale coherent structures like jets and vortices. 
Jovian atmosphere shows a wide diversity of structures: 
Jupiter's great red spot, white ovals, brown barges,... 
One question of fundamental interest is to understand 
and predict the structure and the stability of these quasi 
stationary states (QSS). To that purpose. Miller [l| and 
Robert & Sommeria 0] have proposed a statistical me- 
chanics of the 2D Euler equation (a similar statistical 
theory had been developed earlier by Lynden-BcU 3] to 
describe the violent relaxation of coUisionless stellar sys- 
tems governed by the Vlasov equation; see Q for a de- 
scription of this analogy) . The key idea is to replace the 
deterministic description of the flow u;(r, t) by a proba- 
bilistic description where p{r, a, t) gives the density prob- 
ability of finding the vorticity level lo — a in v a,t time 
t. The observed (coarse-grained) vorticity field is then 
expressed as Tlj{T,t) — J pada. To apply the statistical 
theory, one must first specify the constraints attached 
to the 2D Euler equation. The circulation F = J ujdr 
and the energy E = ^ J TDtpdr will be called robust con- 
straints because they can be expressed in terms of the 
coarse-grained field 10 (the energy of the fluctuations 
can be neglected). These integrals can be calculated 
at any time from the coarse-grained field lo{y, t) and 
they are conserved by the dynamics. By contrast, the 
Casimir invariants If = J f{uj)dr, or equivalently the 



fine-grained moments of the vorticity F^'^-^ = J a;"c?r, 
where lu""- = J pa'^da, will be called fragile constraints 
because they must be expressed in terms of the fine- 
grained vorticity. Indeed, the moments of the coarse- 
grained vorticity F^'^j^ = JtU"(ir are not conserved since 
w" 7^ tU" (part of the coarse-grained moments goes into 
fine-grained fluctuations). Therefore, the moments F//^]^ 
must be calculated from the fine-grained field Lu{r,t) or 
from the initial conditions, i.e. before the vorticity has 
mixed. Since we often do not know the initial conditions 
nor the fine-grained field, the Casimir invariants often 
appear as "hidden constraints" . 

The statistical theory of Miller-Robert-Sommeria 
(MRS) is based on three assumptions: (i) it is assumed 
that the evolution of the flow is strictly described by the 
2D Euler equation (no forcing and no dissipation); (ii) it 
is assumed that we know the initial conditions (or equiv- 
alently the value of all the Casimirs) in detail; (iii) it is 
assumed that mixing is efficient and that the evolution is 
ergodic so that the system will reach, at statistical equi- 
librium, the most probable (most mixed) state. Within 
these assumptions [17], the statistical equilibrium state 
of the 2D Euler equation is obtained by maximizing the 
mixing entropy 



S[p] 



plnp drda, 



(2) 
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at fixed energy E and circulation F (robust con- 
straints) and fixed fine-grained moments ^n>i (fragile 
constraints). We must also account for the normaliza- 
tion condition J pda — 1 . This optimization principle is 
solved by introducing Lagrange multipliers, writing the 
first order variations as 0, : 

5S-(35E-a5T^^an5Tl;3 - I ({r)Spdcrdr = 0. (3) 

n>l 

In the MRS approach, the conservation of all the 
Casimirs has to be taken into account. However, in geo- 
physical situations, the flows are forced and dissipated at 
small scales (due to convection in the jovian atmosphere) 



2 



so that the conservation of the Casimirs is destroyed. 
ElUs-Haven-Turkington |6| have proposed to treat these 
situations by fixing the conjugate variables an>i instead 
of the fragile moments T{^S_'^ (this is essentially a sugges- 
tion that has to be tested in practice). If we view the 
vorticity levels as species of particles, this is similar to 
fixing the chemical potentials instead of the total num- 
ber of particles in each species. Therefore, the idea is 
to treat the fragile constraints canonically, whereas the 
robust constraints are still treated microcanonically. A 
rigorous mathematical formalism has been developed in 
7 1 and a more physical presentation has been given in 
^ . In the EHT approach, the relevant thermodynamical 
potential (grand entropy) is obtained from the mixing 
entropy ([2]) by using a Legendre transform with respect 
to the fragile moments [8j: 



(4) 



n>l 



Expliciting the fine-grained moments, we obtain the rel- 
ative (or grand) entropy 



In 



drda, 



(5) 



where we have defined the prior vorticity distribution 
x{(j) = exp{— ^^^-^ Q;„cr"}. We shall assume that this 
function is imposed by the small-scale forcingso it has 
to be given a priori as an input in the theory 0, Q • 



II. EQUILIBRIUM STATISTICAL MECHANICS 
WITH A PRIOR VORTICITY DISTRIBUTION 

When a prior vorticity distribution is given, the statis- 
tical equilibrium state is obtained by maximizing the rel- 
ative (or grand) entropy S-^ at fixed energy E, circulation 
r and normalization condition J pda = 1 (grand micro- 
canonical ensemble). The conservation of the Casimirs 
has been replaced by the specification of the prior x(cr). 
Writing the first order variations as 63^ — f36E — aST — 
J ({r)6pdadr — 0, we get the Gibbs state 



predicts that the coarse-grained vorticity aJ(r) is a sta- 
tionary solution of the 2D Euler equation and that the 
lJ — relationship is a monotonic function which is in- 
creasing at negative temperatures /3 < and decreasing 
at positive temperatures /3 > 0. We have uj'{ip) — —Puj2- 
We note that the ZJ — tp relationship predicted by the sta- 
tistical theory can take a wide diversity of forms (usually 
non-Boltzmannian) depending on the prior xi'^)- Fur- 
thermore, the coarse-grained distribution ([7]) extremizes 
a generalized entropy in oJ-space of the form ^ : 



S[u 



C{Lj)dr, 



(8) 



at fixed circulation and energy (robust constraints). 
Writing the first order variations as SS — (35E — aST — 0, 
leading to 



C"(cJ) = -Pip - a. 



(9) 



and comparing with Eq. ([7]), we find that C'{x) = 
—F^^{x). Therefore, C is a convex function (C" > 0) 
determined by the prior x(cr) encoding the small-scale 
forcing according to the relation 



C(p) = - / F-Hx)dx 



[(Inx)'] ^{-x)dx. 



We have uj'{4') = —P/C"{uj). Comparing with uj'{ip) 
—j3L02 we find that, at statistical equilibrium 



(10) 



U2 = 1/C"(CJ), 



(11) 



P(r,cr) = 



Z(r 



■x(o')e" 



(6) 



which links the centered variance of the vorticity to the 
coarse-grained vorticity and the generalized entropy. It 
also clearly establishes that C" > 0. On the other hand, 
the equilibrium coarse-grained vorticity (U(r) maximizes 
the generalized entropy ([8|)-(ll0p at fixed circulation and 
energy iff p{r, a) maximizes at fixed E, T (see Ap- 
pendix |^ and jg. [loj). Therefore, the maximization of 
S\uj] at fixed E and T is a necessary and sufficient con- 
dition of EHT thermodynamical stability. 

The preceding relations are also valid in the MRS ap- 
proach except that x(cr) is determined a posteriori from 
the initial conditions by relating the Lagrange multipli- 
ers a„>i to the Casimir constraints ^^yi- In this case of 



with Z = x(o-)e-('''^+")'"dcr. This is the product of 
a universal Boltzmann factor by a non-universal function 
x(cr) fixed by the forcing. The coarse-grained vorticity is 
given by 

f X(cr)cre-('3V'+")^rfcr 

/x(cr)e-('^'^+")'^dCT ^ ^ ^ ^ 

with F{^) — — (lnx)'(*&), where we have defined x{^) ~ 
I-^ X{<^)e^'"^da. It is easy to show that F'($) 

-~UJ2{^) < 0, where uj2 ~ '-^^ ^ 
tered variance of the vorticity. Therefore, F{^) is a de- 
creasing function. Since ZU = f{ip), the statistical theory 



freely evolving flows, the generalized entropy ([8|l- (|T0|) de- 
pends on the initial conditions, while in the case of forced 
flows considered here, it is intrinsically fixed by the prior 
vorticity distribution. On the other hand, a maximum of 
SyJip] at fixed E and T is always a maximum of S[p] at 
fixed E, r and r{'^'^ (more constrained problem). There- 
fore, a maximum of the generalized entropy S^] at fixed 
E and T determines a statistical equilibrium state in the 
MRS viewpoint [lo| . However, the converse is wrong in 
case of "ensemble inequivalence" [ill, [13 with respect to 
the conjugate variables (F^'^-^, a„) because a maximum of 
S[p\ at fixed E, T and F^'^ j^ is not necessarily a maximum 
> is the local cen- of S'j^[p] at fixed E and F. Therefore, the maximization of 



S[uj\ at fixed E and T is a sufficient (but not necessary) 
condition of MRS thermodynamical stability. 
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III. RELAXATION TOWARDS EQUILIBRIUM 

In the case where a smah-scale forcing imposes a prior 
vorticity distribution x(''')i it is possible to propose a 
thermodynamical parametrization of the turbulent flow 
in the form of a relaxation equation that conserves the 
circulation and the energy (robust constraints) and that 
increases the generalized entropy ([8|-(|10|) fixed by the 
prior. This equation can be obtained from a general- 
ized Maximum Entropy Production principle (MEPP) in 
cU-space j9i] . We write lu = to + Cb and take the local aver- 
age of the 2D Euler equation ((!]). This yields DZU/Dt = 
— V-wu = — V- J where D/Dt = d/dt + u-V is the mate- 
rial derivative and J is the turbulent current. Then, we 
determine the optimal current J which maximizes the 
rate of entropy production S = — /C"(ZU)J • VZUdr at 
fixed energy E — J 3 ■ Wipdr — 0, assuming that the en- 
ergy of the fluctuations J^/2lZJ is bounded. According to 
this phenomenological principle, we find that the coarse- 
grained vorticity evolves according to [!,[§]: 



dUJ 

— + u • VcJ = V 

ot 



D 



m 



m = - 



uJ = -Alp, 



(12) 



(13) 



where f3(t) is a Lagrange multiplier enforcing the energy 
constraint i? = at any time. It is shown in Q that 
these equations increase monotonically the entropy {H- 
theorem S' > 0) provided that D > 0. Furthermore, a 
steady state of is linearly dynamically stable iff it 
is a (local) entropy maximum at fixed circulation and 
energy (minima or saddle points of entropy are linearly 
unstable). Therefore, the relaxation equations (|12p - (fT51) 
generically converge towards a (local) entropy maximum 
(if there is no entropy maximum the solutions of the re- 
laxation equations can have a singular behaviour). If 
there exists several local entropy maxima the selection 
will depend on a complicated notion of basin of attrac- 
tion. The diffusion coefficient D is not determined by 
the MEPP but it can be obtained from a Taylor's type 
argument leading to D = Ke^uy^ where e is the coarse- 
graining mesh size and X is a constant of order unity 
[Sy. Assuming that the relation piD remains valid out- 
of-equilibrium (see Appendix C of [8]), we get the closed 
expression D = Ke^ / ^ C" (tU) . This position depen- 
dant diffusion coefficient, related to the strength of the 
fiuctuations, can "freeze" the system in a sub-region of 
space ( "bubble" ) and account for incomplete relaxation 
and lack of ergodicity [1, [l^. The relaxation equation 
(fT2|) belongs to the class of nonlinear mean field Fokker- 
Planck equations introduced in . This relaxation equa- 
tion conserves only the robust constraints (circulation 
and energy) and increases the generalized entropy ([8|)- 
pop fixed by the prior vorticity distribution x((t). It 
differs from the relaxation equations proposed by Robert 



& Sommeria ^li*! for freely evolving flows which conserve 
all the constraints of the 2D Euler equation {E, T and all 
the Casimirs) and monotonically increase the mixing en- 
tropy (21). In Eqs. (112)- US]), the specification of the prior 
x((7) (determined by the small-scale forcing) replaces the 
specification of the Casimirs (determined by the initial 
conditions). However, in both models, the robust con- 
straints E and F are treated microcanonically (i.e. they 
are rigorously conserved). The relaxation equations of 
Robert & Sommeria [l3| and Chavanis Q are essentially 
phenomenological in nature but they can serve as numer- 
ical algorithms to compute maximum entropy states. In 
that context, since we are only interested by the station- 
ary state (not by the dynamics), we can take D — Cst. 
and drop the advective term in the relaxation equation. 
Then, Eq. (|12p can be used to construct (i) arbitrary 
EHT statistical equilibria (ii) a subset of MRS statistical 
equilibria (see the last paragraph of Sec. [IT)) . 



IV. EXPLICIT EXAMPLES 

Let us consider, for illustration, the prior vorticity dis- 
tribution x{^) introduced by EUis-Haven-Turkington Q 
in their model of jovian vortices. It corresponds to a 
de-centered Gamma distribution 



1 



R 



1 



1 



(14) 



where R{z\ a) = F(a) 



-1 .^a-l. 



for z > and R 



otherwise. The scaling of xif^) is chosen such that (a) ~ 
0, var(a) = {a^) = ^2 and skew(a) = {a^)/{a^y^/^ = 
2nl^^X. We get 



(l + Af^2$)i/("^^' 



S7(<I>) = -(lnx)'(<i>) = 



1 + Af^2*' 

Inversing the relation (|16p. we obtain 

1 UJ 



f72 1 + AtU 



(15) 



(16) 



(17) 



After integration, we obtain the generalized entropy 



UJ — — ln(l + Xlu) 
A 



(18) 



In the limit A ^ 0, the prior is the Gaussian distribution 



1 



(19) 



V2^2 

and we have Z($) = e*"^*', ZJ($) = -^^2*, C{uj) = 
The generalized entropy S — J uj'^dr asso- 

ciated with a Gaussian prior is proportional (with the 
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opposite sign) to the coarse-grained enstrophy: S = 
_pc.gy tB]- This gaussian prior leads to Fofonoff 
flows [1^ that have oceanic appHcations. 

When the prior is given by Eq. p^ . the generahzed 
entropy satisfies C"{uj) — l/[ri2(l + ASJ)^] and we obtain 
a parametrization of the form 



duj 
'dt~ 



m 



-u-VtU = V-iD 



VaJ-|-/3(i)f72(l-l-AcU)^VV' 



(20) 



/ DVuj ■ Vipdr 



/L>l]2(l + AcJ)2(VV')2dr 
For A = (Gaussian hmit), we get 

Vlj + /3(i)f72Vi/) 



^ + u • Vw = V • <i D 

ot 



Pit) 



D = Ke^a 



/Dr22(VV')2dr' ^ "—2 



1/2 



(21) 
(22) 
(23) 



Since D and are uniform, we have Duj/Dt = D(Aw — 

P{t)n2uj) with /3(t) = -r2^-(t)/(2f72£:) = ^(t)/^; (to 

arrive at this resuh, we have used integration by parts in 
Eq. 

When the prior has two intense peaks xi'^) — ^{'^ ~ 
"■q) + '5((T — <Ti), the equihbrium coarse-grained vorticity 
is 

^ = '"1+ l + e(-o-.OW+a)- (24) 

This is similar to the Fermi-Dirac statistics. Inverting 
this relation to express <i> = fSif) -I- a as a function of lJ 
and integrating the resulting expression, we obtain the 
generalized entropy 



Sp] = - j [p\np + (1 - p) ln(l - p)]dr. 



(25) 



where uJ = pao + {\ — p)ai. At equilibrium, we have 
L02 ~ l/C"(tU) = (ao — cU) (tU — CTi ) . For the two-peaks 
distribution, we get a parametrization of the form 

+ u- VtU^ V- [D{Vuj + P{t){(Jo-uj){uj - ai)^^})]. 

(26) 



dt 



m = 



J DVlJ ■ V-ipdr 



D = K^JI\ 



(27) 

These are the same equations as in the MRS theory in 
the two levels case w G {f^o,o-i} [11, IS II. They amount 
to maximizing the Fermi-Dirac-like entropy ()25p at fixed 
circulation and energy. This entropy has been used by 
Bouchet & Sommeria [1^ to model jovian vortices. In the 
MRS viewpoint, this entropy describes the free merging 
of a system with two levels of vorticity ctq and a\ while in 
the viewpoint developed here, it describes the evolution 
of a forced system where the forcing has two intense peaks 
described by the prior x(c) = ^(c — co) + ^(c — Ui) 
Q. Other examples of prior vorticity distributions and 
associated generalized entropies are collected in Q. 



V. NONLINEAR DYNAMICAL STABILITY 

Let us consider the Casimir functionals S\u}\ = 
— / C(w)(ir where C is any convex function (C" > 0). 
Since 5, E and F are individually conserved by the 2D 
Euler equation, the maximization problem 



max{S'[w] I E{uj\ = E,Y\uj\ = F} , 



(28) 



determines a steady state of the 2D Euler equation that is 
formally nonlinearly dynamically stable @. Writing the 
first variations as SS — (3SE — aST = 0, the steady state 
is characterized by a monotonic relation lu = F((3ip + 
a) — /(?/;) where F{x) = {C')^^{—x). Let us introduce 
the Legendre transform J = S — f3E and consider the 
maximization problem 



max{ J[w] = S[lu] - l3E[uj] \T[lj] = F} , 



(29) 



If we interpret J as an energy-Casimir functional, the 
maximization problem (|29[) corresponds to the Arnold 
criterion of formal nonlinear dynamical stability. The 
variational problems (|28p and have the same crit- 
ical points (cancelling the first variations) but not nec- 
essarily the same maxima (regarding the second varia- 
tions). A solution of (|29p is always a solution of the 
more constrained problem (|28p . However, the reciprocal 
is wrong. A solution of is not necessarily a solution 
of ([29]) . The maximization problem ((29l) , and the associ- 
ated Arnold theorems, provide just a sufficient condition 
of nonlinear dynamical stability. The criterion of 
Ellis-Haven-Turkington is more refined and allows to con- 
struct a larger class of nonlinearly stable steady states. 
For example, important equilibrium states in the weather 
layer of Jupiter are nonlinearly dynamically stable ac- 
cording to the refined stability criterion while they 
do not satisfy the Arnold theorems @. The maximiza- 
tion problem (|29p determines a subclass of solutions of 
the maximization problem (|28p . This is similar to a sit- 
uation of "ensemble inequivalence" with respect to the 
conjugate variables {E,f3) in thermodynamics [HI, [l^ . 
Indeed, (|28p is similar to a criterion of "microcanonical 
stability" while ([^5]) is similar to a criterion of "canonical 
stability" in thermodynamics, where S is similar to an en- 
tropy and J is similar to a free energy Q . Canonical sta- 
bility implies microcanonical stability but the converse 
is wrong in case of ensemble inequivalence fl8'|. Since 
the relaxation equations (|12p - (|13p solve the maximization 
problem (|28p . they can serve as numerical algorithms to 
compute nonlinearly dynamically stable stationary solu- 
tions of the 2D Euler equation according to the criterion 
of Ellis-Haven-Turkington. Note that if we fix /3, the re- 
laxation equation increases monotonically the "free 
energy" J — S — I3E (i/-theorem, J > 0) until a (local) 
maximum of J at fixed F is reached [^. Therefore, we ob- 
tain a numerical algorithm that solves the maximization 
problem (|29p and determines a subclass of nonlinearly 
dynamically stable stationary solutions of the 2D Euler 
equation corresponding to the Arnold criterion. 
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VI. CONCLUSION 

In this paper, we have shown that the maximization 
of a functional S[u!] at fixed circulation T and energy E 
in 2D turbulence can have several interpretations. When 
S is given by ([8l)- (|10p . this maximization problem de- 
termines: (i) The whole class of stable EHT statistical 
equilibria for a given prior vorticity distribution xi'^) 
fixed by the small-scale forcing, (ii) A subclass of stable 
MRS statistical equilibria for initial conditions leading 
to a vorticity distribution xi'^) statistical equilibrium. 
When S is given by ^ where C is an arbitrary convex 
function, this maximization problem determines a non- 
linearly dynamically stable stationary solution of the 2D 
Euler equation according to the refined EHT criterion. 
The next step is to determine whether particular forms 
of generalized entropies are better adapted than others 
to describe specific flows and whether they can be re- 
grouped in "classes of equivalence" [9j]. For example, the 
enstrophy functional turns out to be relevant for certain 
oceanic situations [l^ and the Fermi-Dirac-like entropy 
for jovian flows [11]. Working with a suitable generalized 
entropy S[uj] with only two constraints (T, i?) is more 
convenient than working with an infinite set of Casimirs 
as in the MRS theory. This reduced maximization prob- 
lem is still very rich because, for any considered form 
of generalized entropy S[uj], many bifurcations can take 
place in the parameter space {E,r) [sl. [l6j . 

APPENDIX A: GENERALIZED ENTROPY 

We can introduce the generalized entropy S\lJ] in 
the following manner. Initially, we want to determine 
the vorticity distribution /9*(r, a) which maximizes S^[p\ 
with the robust constraints E\lJ] = E, r[ZU] = F, and the 
normalization condition J pda ^ 1. To solve this max- 



imization problem, we can proceed in two steps. First 
step: we determine the distribution pi(r, cr) which max- 
imizes [p] with the constraints J pda = 1 and a fixed 
vorticity profile J pa da ~ uj{r) (note that fixing ZJ auto- 
matically determines F and E) . This gives a distribution 
pi[uj(r), a] depending on Lu{r) and a. Substituting this 
distribution in the functional S';(;[p], we obtain a func- 
tional S\lJ] = S^[pi] of the vorticity ZJ. Second step: 
we determine the vorticity field w*(r) which maximizes 
S[ID] with the constraints E[uj] = E and T[uj] = F. Fi- 
nally, we have p^,{r,a) — pi[iU*(r), cr]. Let us be more 
explicit. The distribution yOi(r, ct) that extremizcs S-^[p\ 
with the constraints J pda = 1 and J pa da — uj{r) satis- 
fies the first order variations SS^ — J ^{r)S{J pada)dv — 
J ({r)5{J pda)dr = 0, where $(r) and C(r) are Lagrange 
multipliers. This yields 

Pi(r,a) = ^x(^)e-^*('-\ (Al) 

where Z{r) and $(r) are determined by 
Z{r) = /x(cr)e-'"*('")dCT = and w(r) = 

^/x(cr)CTe~'"*Wdcr = -(lnx)'($)- This critical 
point is a maximum of with the above-mentioned 
constraints since 6^S^ = — / ^^^dvda < 0. Then 

'S'xIpi] = J pi{a^ + \nx) drda = J {uJ^ + \n x{^)) dr . 
Therefore S\uj] = S^[pi] is given by S\uj] — — J C{lJ) dr 
with C{lJ) = -aJ$ - \nx{^). Now, $(r) is related 
to a;(r) by aJ(r) = — (Inx)'(^'). This implies that 
C"(cJ) = -$ = -[{lnxY]-H-^) so that 

CiU) = - [ilnx)r\-^)dx. (A2) 

This is precisely the generalized entropy (fTO|) . Therefore, 
p» (r, a) = pi [ZZJ* (r) , cr] is a maximum of [p] at fixed E 
and F iff ZU*(r) is a maximum of S\lu] at fixed E and F. 
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